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Abstract We provide lower error bounds for randomized algorithms that approx- 
imate integrals of functions depending on an unrestricted or even infinite number 
of variables. More precisely, we consider the infinite-dimensional integration prob- 
lem on weighted Hilbert spaces with an underlying anchored decomposition and 
arbitrary weights. We focus on randomized algorithms and the randomized worst 
case error. We study two cost models for function evaluation which depend on the 
number of active variables of the chosen sample points. Multilevel algorithms be- 
have very well with respect to the first cost model, while changing dimension algo- 
rithms and also dimension-wise quadrature methods, which are based on a similar 
idea, can take advantage of the more generous second cost model. We prove the 
first non-trivial lower error bounds for randomized algorithms in these cost mod- 
els and demonstrate their quality in the case of product weights. In particular, we 
show that the randomized changing dimension algorithms provided in [L. Plaskota, 
G. W. Wasilkowski, J. Complexity 27 (2011), 505-518] achieve convergence rates 
arbitrarily close to the optimal convergence rate. 



1 Introduction 

Integrals over functions with an unbounded or infinite number of variables are im- 
portant in physics, quantum chemistry or in quantitative finance, see, e.g., (8] |25] 
and the references therein. In the last few years a large amount of research was 
dedicated to design new algorithms as, e.g., multilevel and changing dimension al- 
gorithms or dimension-wise quadrature methods, to approximate such integrals ef- 
ficiently. Multilevel algorithms were introduced by Heinrich in HI 21 in the context 
of integral equations and by Giles in [8] in the context of stochastic differential 
equations. Changing dimension algorithms were introduced by Kuo et al. in lfl6l 
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in the context of infinite-dimensional integration on weighted Hilbert spaces and 
dimension-wise quadrature methods were introduced by Griebel and Holtz in ifTTl 
for multivariate integration; changing dimension and dimension-wise quadrature al- 
gorithms are based on a similar idea. 

Here we want to study the complexity of numerical integration on a weighted 
Hilbert space of functions with infinitely many variables as it has been done in 

iBimnsimiiEDiainm. Hie Hilbert s P ace we c ° nsider here aii ° ws f ° r a 

so-called anchored function space decomposition. For a motivation of this specific 
function space setting and connections to problems in stochastics and mathematical 
finance see, e.g., |[T3l[T8l . We derive lower error bounds for randomized algorithms 
to solve the infinite-dimensional integration problem. Notice that the complexity 
of integration problems is less well understood in the randomized setting than in 
the deterministic setting (where only deterministic algorithms are permitted and the 
deterministic worst case error is considered), see, e.g., the comments in l20l p. 487]. 

Our error bounds are for the randomized worst case error and are expressed in 
terms of the cost of a randomized algorithm. Here we solely take account of func- 
tion evaluations, i.e., the cost of function sampling, and disregard other cost as, e.g., 
combinatorial cost. Notice that this makes the statements of our lower bounds only 
stronger. To evaluate the cost of sampling, we consider two sampling models: the 
nested subspace sampling model (introduced in J5), where it was called variable sub- 
space sampling model) and the unrestricted subspace sampling model (introduced in 
lfl6ll ). Our lower error bounds are the first non-trivial lower bounds in these settings, 
cf. also the comments in the introductions of lfl3l |2T1 . Due to space restrictions, 
we do not provide new constructive upper error bounds. For the same reason we 
refer for a formal definition of multilevel algorithms and changing dimension algo- 
rithms for the infinite-dimensional integration problem on weighted Hilbert spaces 
to fPH [THl l9l and |[T6ll2D . respectively. In this article we only compare our lower 
bounds to already known upper bounds. In particular, we show that the random- 
ized changing dimension algorithms provided for product weights in ||2T1 achieve 
convergence rates arbitrarily close to the optimal rate of convergence. 

Let us mention that similar general lower error bounds for infinite-dimensional 
integration on weighted Hilbert spaces are provided in [6i] in the determistic setting 
for the anchored decomposition and in [4] in the randomized setting for underlying 
ANOVA-type decompositions (to treat the latter decompositions, a technically more 
involved analysis is necessary). 

The article is organized as follows: In Section|2]the setting we want to study is in- 
troduced. In Section|3]we prove new lower bounds for the complexity of randomized 
algorithms for solving the infinite-dimensional integration problem on weighted 
Hilbert spaces. In Section lXTl we provide the most general form of our lower bounds 
which is valid for arbitrary weights. In Section [572] we state the simplified form of 
our lower bounds for specific classes of weights. In particular, we show in Section 
13.2. ll that the randomized changing dimension algorithms from l2TI are essentially 
optimal. 
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In this section we describe the precise setting we want to study. A comparison with 
the (slightly different) settings described in the papers j9][l6l|2T] will be provided in 
the forthcoming paper ifTol ; we refer to the same paper and to lfl4l[T3l for rigorous 
proofs of the statements on the function spaces we consider here. 



2.1 Notation 

For neNwe denote the set {l,...,n} by [n]. If u is a finite set, then its size is 
denoted by \u\. We put ^:={wCN||m|<°°].. We use the common Landau symbol 
O, and additionally for non-negative functions f,g : [0,°°) — > [0,°°) the notation 



2.2 The function spaces 

As spaces of integrands of infinitely many variables, we consider reproducing kernel 
Hilbert spaces; our standard reference for these spaces is 0"). 

We start with univariate functions. Let D C R be a Borel measurable set of R 
and let Ic : D x D — > R be a measurable reproducing kernel with anchor a € D, i.e., 
Ic(a,a) = 0. This implies k(-,a) = 0. We assume that k is non-trivial, i.e., k ^ 0. We 
denote the reproducing kernel Hilbert space with kernel k by H = H{k) and its scalar 
product and norm by (•,•}/? and || • ||#, respectively. Additionally, we denote its 
norm unit ball by B(k). We use corresponding notation for other reproducing kernel 
Hilbert spaces. If g is a constant function in H, then the reproducing property implies 
g = g(a) = (g,k(-,a))H = 0. Let p be a probability measure on D. We assume that 



where by convention k® = 1. The Hilbert space with reproducing kernel k u will 
be denoted by H u = H{k u ). Its functions depend only on the coordinates jeu. 
If it is convenient for us, we identify H u with the space of functions defined on 
D" determined by the kernel Y[jeuk(xj,yj), and write / M (x M ) instead of /«(x) for 
/„ £ H u , x E D N , and x„ := (xj)j eu £ D". For all /„ e H u and x G D N we have 



f = Q(g) if g = 0(f). 




(1) 



For arbitrary x, y € D and ugf we define 



/„ (x) =0 if xj = a for some j <E u. 



(2) 
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This property yields an anchored decomposition of functions, see, e.g., (T7J- 

Let now y = (y^ue'W be weights, i.e., a family of non-negative numbers. We 
assume 

£ y u M H < 00. (3) 
Let us define the domain 2£ of functions of infinitely many variables by 



E Yuku(xj,Xj)< 00 

«c -U 1 



Let be the infinite product probability measure of p on Z) N . Due to our assump- 
tions we have = 1, see |[T3l Lemma 1] or ifTUl . We define 

K y (x,y):= E Yukufaj) forallx,y G 9£ . 

Ky is well-defined and, since Ky is symmetric and positive semi-definite, it is a 
reproducing kernel on SC x S£ , see fl~|- We denote the corresponding reproducing 
kernel Hilbert space by My = H(Ky), its norm by || • |L, and its norm unit ball by 
By = B(Ky). For the next lemma see 03 Cor. 5] or IfTUl . 

Lemma 1. The space My consists of all functions f = Y.ue"?/ fu> fu G H u , that have 
a finite norm 

= (Lru 1 \\Mk) ■ 



7- 



For u G let P u denote the orthogonal projection P u : My —> H u , f H > f, onto 
H u . Then each / G My has a unique representation 

/= E /„ with /„ = P u (f) G //„, U G flr. 



2.3 Infinite-dimensional integration 

For a given / G My we want to approximate the integral 



/(/):= y £ ./(x)M(dx). 
Due to (01, / is continuous on My and its representer h G J^y is given by 



h(x) = J^K r {x,y)n(dy). 
The operator norm of the integration functional / is given by 



Lower Error Bounds for Randomized Multilevel and Changing Dimension Algorithms 5 

\\n 2 j? Y =n 2 r = e y«c 1 o 1 <°°, (4) 

ue°7/ 

where 

C := / / *(jc,y)p(dje)p(dy). 
Jd Jd 

We have Co < M. We assume that / is non-trivial, i.e., that Co > and y u > for at 
least one «ef . For uGf and / G we define /„ .= IoP u , i.e., 



/«(/)= / / H (x H )p"(d*.«), 
and the representer ft M of /„ in M'y is given by h„(x) = P u (h)(x). Thus we have 
h=Y,h u and /(/)=£/„(/„) for all /GJ^. 

Furthermore, 

\\hu\\l=JuC^ for all (5) 



2.4 Randomized algorithms, cost models, and errors 

As in 1131 . we assume that algorithms for approximation of /(/) have access to the 
function / via a subroutine ("oracle") that provides values /(x) for points x G D N . 
For convenience we define /(x) = for x G D N \ 3£ ' . 

We now present the cost models introduced in @ and fl6l : In both models we 
only consider the cost of function evaluations. To define the cost of a function evalu- 
ation, we fix a monotone increasing function $(v) : No — > [1,°°]- For our lower error 
bounds we will later assume that $( v) = Q (V s ) for some s > 0. For each vG^we 
define the finite-dimensional affine subspace 3£ v>a of 3£ by 

X V)ll := {x G Z) N \xj = a for all j G N\v}. 

In the nested subspace sampling model function evaluations can be done in a 
sequence of affine subspaces 5£ Vl A C 5£ V2 >a C ■ ■ ■ for a strictly increasing sequence 
v = (vj);-gN °f sets 7^ G % , and the cost for each function evaluation is given by 
the cost function 

c V)a (x):=inf{$(|vi|)|xG^r Vl , fl }, (6) 

where we use the standard convention that inf0 = °°. Let C ne st denote the set of all 
cost functions of the form ©. The nested subspace sampling model was introduced 
in 0, where it was actually called "variable subspace sampling model". We prefer 
the name "nested subspace sampling model" to clearly distinguish this model from 
the following cost model, which is even more "variable": 
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In the unrestricted subspace sampling model we are allowed to sample in any 
subspace 2£ u ^ a , «6f , without any restriction. The cost for each function evaluation 
is given by the cost function 

e (x):=inf{$(|H|)|xe < r H , fl ,iie < ar}. (7) 

Put C unr := {c a }. The unrestricted subspace sampling model was introduced in |[T6ll . 
where it did not get a particular name. Obviously, the unrestricted subspace sam- 
pling model is more generous than the nested subspace sampling model. 

We consider randomized algorithms for integration of functions / £ Jtfy. For a 
formal definition we refer to (5] |T9] [23] [24] . Here we require that a randomized al- 
gorithm Q yields for each / £ Jtfy a square-integrable random variable Q(f). (More 
precisely, a randomized algorithm Q is a map Q : Q x Jtfy — >• R, (ft),/) n- Q(co,f), 
where £2 is some suitable probability space. But for convenience we will not specify 
the underlying probability space £2 and suppress any reference to Q or co £ Q . We 
use this convention also for other random variables.) Furthermore, we require that 
the cost of a randomized algorithm Q, which is defined to be the sum of the cost 
of all function evaluations, is a random variable, which may depend on the function 
/. That is why we denote this random variable by cost c (g,/), c the relevant cost 
function from C nes t or C unr . 

We denote the class of all randomized algorithms for numerical integration on 
Jtfy that satisfy the very mild requirements stated above by stf™. For unrestricted 
subspace sampling we additionally consider a subclass g/ les of £/ lan . We say that an 
algorithm Q £ si/ lan is in £/ les if there exist an n £ No and sets vi, . . . , v„ £ ^ such 
that for every / £ My the algorithms Q performs exactly n function evaluations of /, 
where the zth sample point is taken from 3€v ua , and E(cost Co (Q,f)) = L/'=i $(l v i'D- If 
additionally |vi|, . . . , \v n \ < (0 for some CO £ N, we say that Q £ £/ Tes ~ m . Notice that 
the classes stf™, £/ les , and £/ les ~ m contain in particular non-linear and adaptive 
algorithms. 

The worst case cost of a randomized algorithm Q on a class of integrands F is 

cost nes t(2,^) := inf supE(cost c (2,/)) 

in the nested subspace sampling model and 

costunrCG.F) := supE(cost, fl (Q,f)) 

feF 

in the unrestricted subspace sampling model. The randomized (worst case) error 
e(Q,F) of approximating the integration functional / by Q on F is defined as 

e(Q,F):= (jupE ((/(/)- Q(f)) 2 

For JVgR, mod £ {nest,unr}, and * £ {ran, res, res — co} let us define the corre- 
sponding Nth minimal error by 
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2.5 Strong Polynomial tractability 

Let oeN, mod £ {nest,unr}, and * G {ran, res, res — 03}. The e-complexity of the 
infinite-dimensional integration problem / on Jtify in the considered cost model with 
respect to the class of admissable randomized algorithms &/* is the minimal cost 
among all admissable algorithms, whose randomized errors are at most e, i.e., 

comp^ od ( £ ,B r ) := inf {cost mod (<2,B r ) | Q E stf* and e(Q,B Y ) < e} . (8) 

The integration problem / is said to be strongly polynomially tractable if there are 
non-negative constants C and p such that 

comp^ od (e,B r ) < Ce~ p forall£>0. (9) 

The exponent of strong polynomial tractability is given by 

PmoA = Pmod(y) : = lv ^{p I P satisfies for some C > 0}. 

Essentially, 1 / p* mQi is the convergence rate of the Mh minimal error e^ od (A^,By). In 
particular, we have for all p > p* moA that e^ od (^V,B r ) = 0(^V- 1 /"). 



3 Lower Bounds 

We start in Section lTTI bv proving lower bounds for general weights. In Section [3T2l 
we show how these bounds simplify for several specific classes of weights. 



3.1 Results for General Weights 

Let 7 = (YuJue'W be a given family of weights that satisfy ([3j. We denote by y the 
family of weights defined by 

y«:=y«C^ for all sef. (10) 

Recall that (O implies £ u6 ^ % < 00 . Weights y are called finite-order weights of 
order co if there exists an(i)£N such that y, = for all u 6 % with |k| > 03, Finite- 
order weights were introduced in Q for spaces of functions with a finite number of 
variables. The following definition is taken from (9|. 

Definition 1. For weights y and c £ N let us define the cut-off weights of order a 
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I otherwise. 

Cut-off weights of order a are in particular finite-order weights of order a. Let 
us denote by u\(o),U2{o), ■ ■ ■, the distinct non-empty sets «6f with yi' 7 ' > for 
which Y^) a \ > h£(a) — ^ et us P ut M o(c) := 0- We can make the same definitions 
for a = °°; then we have obviously y^°°) = y. For convenience we will often suppress 
any reference to a in the case where a = °o. For (J6NU {°°} let us define 



decay y a := sup jp G R | lim y^/ = 



Due to assumption (01 the weights we consider always satisfy decay ^ a > 1 for all 
(JSNU {°°}. The following definition is from |9l . 

Definition 2. For a e NU {°°} let t* e [0,«s] be defined as 

t* :=M{t >0| 3Q >0Vv€ f/ : |{/ £ N|m,-(ct) C v}| <C,|v| r }. 

Let a e N. Since |m,-((t)| < o for all i e N, we have obviously f* < a. On the 
other hand, if we have an infinite sequence (kj(ct)) j e f$, it is not hard to verify that 
t* a > 1, see @. 

For vi, . . . , v n € & we use the short hand {v,} for (v,)" =1 . Put v := U" =1 v,- and 
define the mapping 

;';3('6[n]:H^Cvi 

^{vi},a i s me orthogonal projection of onto the subspace 

H {vi},a : = L H "j- 

j;3ie[n]:ujCvi 

Put 

b {l , } ,„:= sup |7(/)-/(f {v . }ia /)|. 
feBy 

In the case where n = 1 and v = Vi, we simply write % a and b vfl . In that case we 
have, due to (0, 

(S / V)a (/))(x)=/(x v ;a) for all / e Jf r and x e JT, (13) 
where the jth component of (x,,;a) is defined by 



(x v ;a); := 

Lemma 2. Let Vi , . . . , v n £ Then 



xj ify'ev, 

a otherwise. 
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Proof. Let fi! v .y a denote the representer of the continuous functional / o f{ Vl j a . Due 
to (O we get 

h {M,« = L h "j- 

Since /j — h{ Vi },a is the representer of / — / o f{ Vi } a in j?^, we obtain with (f5]) 

2 



b{v,-}, fl = l|/'-/Hv,}, a llr 



j-yie[n]:Uj<^Vi 



j;Vie[«]:«j^v,- 



y;Vie[n]:«;^vj 



Lemma 3. Lef G (1 /2, 1] and vi, . . . , v n G Lef f/ie randomized algorithm Q G 
^ ran M f/s/> P(g(/) = C2(«P {vj}ifl /)) > 9 for all f G B r . Then 

e{Q,B 7 f>{2e-\)b\ vtU . 

Proof. Since fr,,.} a is an orthogonal projection, we have for all / G By that g := 
f - V{ Vi },af G By. Furthermore, f {v , } , a (s) = f{ v ,}, fl (-§) = 0. Let A := - 
Q(-g)}. Then {2(g) = e(f {v . } , a g)} D {Q(-g) = Q(V {vi} , a (-g))} C A, and hence 
P(A) > 20- 1. Therefore 

e(Q,B 7 ) 2 >max {e ((/(g) - Q(g)) 2 ) ,E ((/(-*) - G(-^))) 2 ) } 

/ (I(g)-Q(g)) 2 dP,[ (I(-g)-Q(-g))) 2 d 

J A J A 

>(20 - l)|/(g)| 2 = (20 - l)|/(/) -/(f { ,, v ,/)| 2 . 

Hence e{Q,B y ) 2 > (20 - l)sup /eBy |/(/) -I(W {v . u f)\ 2 = (20 - l)b| v . }|fl . 

Further Assumptions. We assume for the rest of this article that $(v) = £2 (v s ) for 
some s G (0,°°). Furthermore, we assume that y^j > and that there exists an a > 
such that for univariate integration in H(y^k) the Mh minimal error satisfies 



>max 



e IW (N,B(r {l} k)) = Q(N- a / 2 ). 



(14) 



(Note that in the univariate situation the nested and the unrestricted subspace sam- 
pling models are equal; that is why we suppress the reference to unr or nest.) Since 
(/{!}£) G By, assumption ([Pil l implies in particular 

e™(N,B r ) = £2(N- a f 2 ) and ^ a (N,B r ) > e™(N,B y ) = f2(AT«/ 2 ). (15) 

Theorem 1. Consider the nested subspace sampling model. To achieve strong poly- 
nomial tractability for the class jz/ ran it is necessary that the weights satisfy 
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decay y a > 1 for all ffeN. (16) 

If ( 1761 ) holds, we have 

/?™>maxi ^,sup- 



x « aeN decay y a - 1 J 

As we will see in Section [3721 for product weights and finite-order weights con- 
dition ( fT6l ) is equivalent to decay y = decay yoo > 1. 

Proof. Let <2 G j^ ran with cost nest (2,By) < N. Then there exists an increasing se- 
quence v = (v,), GN , 7^ v ( - G ^T, such that E(cost Cvil (g,/)) </V + 1 for every f £B y . 
Let ra be the largest integer satisfying $(|v m |) < 4(/V+ 1). This implies for all / G By 
thatP(2(/) = Q{% m , a f)) > 3/4, see <|T3}. Due to Lemma E] and [3] we get 

e{Q,B y ) 2 > 1 - £ ?„.. 

Let us now assume that 7 are weights of finite order ft). Then we get for t > and 
a suitable constant Q > 

T m := |{;>> C v m }| < Q\v m \ 1 = 0{N'l s ), 

since N = Q(\v m \ s ). Hence we get for p m > decay y a = decay y > 1 

j=x,„+i 

For general weights 7, a G N, and cut-off weights y'' 7 ' we have e{Q,By) > e(Q,B^„)), 

see also J5] Remark 3.3]. Since the cut-off weights yW are weights of finite order 
<7, we get for all /;><j > decay y a and f CT > f * 



e (e,B r ) z =i2(Ar- ( ™J. (17) 

Since (fl3T l holds, the inequality for the exponent of tractability follows. 

Now assume that the infinite-dimensional integration problem I is strongly 
polynomially tractable. Let a G N. Then we get from inequality $1% that p a > 
l+2V(fa// n a e n st ). Hence 

decay y , a >l + -^. 

"nest 

Thus we have decay y ^ > 1 for all a G N. 

Theorem 2. Consider the unrestricted subspace sampling model. To achieve strong 
polynomial tractability for the class £/ ves it is necessary that the weights satisfy 

decay y a > 1 for all a G N. 
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If this is the case, we have 

p UM .>ma X |-,sup^-— — r |. 

Proof. Let Q G £/ les have cost um (Q,By) < N. Then there exists anngN and coordi- 
nate sets vi, . . . , v n such that Q selects randomly n sample points xi G 3%i lt a, ■ ■ ■ ,x„ G 
&v n ,a and £? =1 $(| v;|) < N. Since = Q(f{v r },a/) for all / G B y , we obtain from 
Lemma|2]and|3] 

e(Q,B r ) 2 > £ y„ r 

;;Vie[n]:«^''i 

Let us first assume that y are weights of finite order co. Then we get with Jensen's 
inequality for t > and suitable constants C ( , c > 

|{j I 3i G [n] : C v,}| < £ |{; | M , C v,}| < £Q|v,|' 
i=i i=i 

/ „ \ l/min{l..v/f} 

< Qy L I v,- I s j < c t ( cN ) 1 7 min{ 1 ' s/t} ■ 

Hence we obtain for S := [C,(cA0 1/min{1 '' !/ ' } l and all p ra > decay y a) 

e(g,B r ) 2 > £ =i2(5 1_/ ' ffl ) =i2 (a^s 1 ^ ) . 

i=S+i V / 

If we have general weights y, then we obtain for a G N and the cut-off weights y^ a ' 
that e(Q,B(Ky)) > e(Q,B(Kj a ))). From this and ([15]) the inequality for the expo- 
nent of tractability follows. Similarly as in the proof of Theorem[T] the necessity of 
condition ( TToT ) is easily established. 

Theorem 3. Let CO G N be fixed. We have for the exponent of tractability p^~ m in 
the unrestricted subspace sampling setting 

\ 2 1 
Pum a >max<^ -,sup 



^ a a<m decay yCT - 1 J 

Proof. We follow the lines of the proof of Theorem [2] and use the same nota- 
tion. The difference is that this time Q selects randomly « sample points xi G 
5£ Vl>a , . . . ,x„ G 3£ Vn ,a, where |v,-| < CO for all ; G [«], and that we therefore can make 
the estimate |{; | 3i G [n] : uj C v;}| < 2 m n = 0(N), since $(|v,|) > 1 for all i G [n] 
by definition of the function $. Hence we get this time for p > decay ^ 

e(Q,B r ) 2 > £ y Uj =Q{N l -P). 

j=2 a n+\ 
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A comparison of Theorem [2] and [3] indicates that there are cost functions and 
classes of finite-order weights for which changing dimension algorithms cannot 
achieve convergence rates that are arbitrarily close to the optimal rate. Let us re- 
call that for weights of finite order co, changing dimension algorithms as defined 
in lfl6l Proof of Thm. 5] would only use sample points from sample spaces 3£ ufi 
with \u\ < co; see also the comment at the beginning of Section 4 in lfl6l . Exam- 
ples of such cost functions and finite-order weights would be $(£) = £2(k s ) and 
lexicographically-ordered weights of order co > s, see Section [3.2.3l (A similar ob- 
servation was made for the deterministic setting, see ||9] Thm. 3.2 & Sect. 3.2.3 ].) 



3.2 Results for Specific Classes of Weights 

Here we consider some example classes of weights and show how our bounds from 
Section [3~T1 simplifv in those settings. 

3.2.1 Product weights and finite-product weights 

Definition 3. Let (Yj)jen be a sequence of non-negative real numbers satisfying 

Yi > 72 > With the help of this sequence we define for CO G N U {°°} weights 

Y = (Yu)uev by 



where we use the convention that the empty product is 1 . In the case where co = °°, 
we call such weights product weights, in the case where CO is finite, we call them 
finite-product weights of order CO. 

Product weights were introduced in |f22| and have been studied extensively since 
then. Finite-product weights were considered in J9). Observe that for a G N the 
cut-off weights of product weights y are finite-product weights of order a. 

Let us assume that y are product or finite-product weights. As shown in J5] 
Lemma 3.8], we have 



(Actually, J5] Lemma 3.8] states inequality (fl~9b only for all a G N. But the proof 
provided in |9| is also valid for the case a = °°.) In particular, we see that for strong 
polynomial tractability with respect to the nested subspace sampling model and the 
class £/ lm or with respect to the unrestricted subspace sampling model and the 
class <2/ res it is necessary that decay^, = decay yoo > 1. Since f* = 1, we obtain from 
Theorem[Q|2] and|3] 




decay r j = decay r a for all a G N U {°°}. 



(19) 
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^ st > m ax{|,^— T }, (20) 

and 

res I 2 2min{l,s} 1 m I 2 2 

^ maX {a'd^a^j' ^ ^ maX {«'^— I)' (21) 

Note that the bounds for finite-product weights are the same as for product weights. 

Remark 1. For product and and order-dependent (POD) weights {y u )ue' J 7/^ which 
were recently introduced in lfl"5l and are of the form 

Yu = r ]u \Y[Yj, where 7i > Ji > • • • > 0, and T = H = 1, r z ,r 3 ,... > 0, 

jeu 

identity ( fT9b still holds; for a proof see J6|. Thus (l20l and (fJTJ are also valid for 
POD weights. Product and finite-product weights are, in particular, POD weights. 

Let us assume that there exist constants c,C, K > 0, (X\ > 0, and £ [0, 1] such 
that for all 6 ^ u 6 ^ and all n > 1 there exist randomized algorithms Q n ll using 
for all /„ e H u at most n function values of /„ with 

rCW / \n(n + U \ "id"!" 1 )" 2 

HUfu)-Q,M)\ 2 ) < [ 1 + (J^ J T)k) 

Note that necessarily K < a. Let us further assume that decay y j > 1 and the cost 
function $ satisfies $(v) = 0(e rV ) for some r > 0. Plaskota and Wasilkowski proved 
in lETI with the help randomized changing dimension algorithms that 

^ S ^ maX {l decay 2 yJ -l }- 

Hence, if £2 (v) = $(v) = 0(e rv ) and K = a, our lower bound (f2Tb is sharp and the 
randomized algorithms from H21H exhibit essentially the optimal convergence rate. 

Let us consider a more specific example, namely the case where D = [0, 1], k 
is the Wiener kernel given by k(x,y) = min{x 7 y}, and p is the restriction of the 
Lebesgue measure to D. In this case the anchor a is zero. The space H(k) is the 
Sobolev space anchored at zero, and its elements are the absolutely continuous 
functions / with /(0) = and square-integrable first weak derivative. It is known 
that K = 3 = a, see EU Ex. 1 and Prop. 3] (or ED Ex. 2]) and HU Sect. 2.2.9, 
Prop. 1]. Thus the upper bound from 12T1 and our lower bound OTT i establish for 
Q(v) = $(v) = 0(e rv ) that 

^ = maX {lde^l}- 
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For the same specific example Hickernell et al. showed for the case $(v) = @(v) 
with the help of multilevel Monte Carlo algorithms that 

pS ^ maX { 2 ' decay 2 y[1 -l } for > ^ 
see JT3] Cor. 5]. Hence 

pa = s -^_ T for decay,, G (1,2], 

Similarly as in the deterministic setting |[T8l |6) or in the randomized setting with 
underlying ANOVA-type decomposition JfQ, our lower bound for p™ t is sharp for 
sufficiently large decay y j . This may be proved by using multilevel algorithms based 
on the integration algorithms provided in fl26l Sect. 4] (cf. also lETl Sect. 3.2]) or 
on scrambled scrambled quasi-Monte Carlo algorithms similar to those discussed in 
|3l , but providing a rigorous proof for this claim is beyond the scope of this article. 

3.2.2 Finite-Intersection Weights 

We restate Definition 3.5 from J9)- 

Definition 4. Let p G N. Finite-order weights (y«) M <E^ are called finite-intersection 
weights with intersection degree at most p if we have 

|{ve^|y v >O,«nv^0}| < 1+p for all u G W with Ju > 0. (22) 

For finite-intersection weights of order CO it was observed in J9J that t * = 1 for 
all (J G N, resulting in the lower bounds < [2Qb and (f2TT > with decay y l replaced by 
decay yfi) . 

3.2.3 Lexicographically-ordered weights 

To every set u C N with \u\ = I we may assign a word (p(u) := i\i<i . . . i(, where for 
j G [£) the number ij is the y'th-largest element of u. On the set of all finite words over 
the alphabet N we have the natural lexicographical order -<i ex , where by convention 
the empty word should be the first (or "smallest") word. 

Definition 5. We call weights y lexicographically-ordered weights of order 03 if 
70 = L J u > for all u C N with \u\ < CO, and 

<P(ui) -<iex <p(uj) for a H hi € N satisfying ; < j. 

Lexicographically-ordered weights were introduced in |9 l|. Their properties com- 
plement the properties of the other classes of weights considered before, see ||6] [9] 
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for more information. For lexicographically-ordered weights of order co we have 
f£ = min{(7, co}. Hence we get from Theorem[T||2] and[3]the lower bounds 

D - >maJ- 25/10 \ 
Pnest - maX \a'decay 7i£0 -l/' 

and 

P™^{a> decay y , ffl -l /' ^ " j a ' decay^ - ! J ' 

The lower bounds indicate that in the setting where CO > s and decay y m is only mod- 
erate, changing dimension algorithms (which are algorithms of the class g/ les ~ m ) 
cannot achieve the optimal rate of convergence and can be outperformed by multi- 
level algorithms (which can exhibit a behavior similar to the lower bound for p™ t 
above). For the deterministic setting and the Wiener kernel k(x,y) = min{x,y} on 
D = [0, 1] this was rigorously proved in j9| by lower bounds for changing dimen- 
sion algorithms and upper bounds for multilevel algorithms, see J9] Thm. 3.2 & 
Sect. 3.2.3 ]. 
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